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o : Abstract 

We describe how the complete solution to the two-dimensional constant quantum 
Yang-Baxter equation was found [1-3]. 
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In this talk we give a brief description of how the two-dimensional constant quantum Yang- 
Baxter equation (YBE) was solved. The results were announced in [1], a detailed account 
of the solution process will be published elsewhere [2], see also [3]. 
First let us recall the spectral parameter dependent form of YBE: 
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Here summation over the repeated k indices is understood. The constant quantum YBE 
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is obtained from (1) with u = v = Qoiu = v = ±00. 

The YBE appears in many physical contexts. It was first derived in the study of solv- 
able vertex models in statistical mechanics as the condition of commuting transfer matrices 
(Yang, Baxter), another derivation follows from the factorization of the S-matrix in 1+1 di- 
mensional Quantum Field Theory (Zamolodchikov). The YBE is also essential in Quantum 
Inverse Scattering Method for integrable systems as developed by the Leningrad school. The 
more abstract setting is in terms of Sklyanin algebras, quantum groups and Hopf algebras 
(Drinfeld). More recently the connection to braid groups and knot theory has been stud- 
ied, here the spectral parameter independent form arises naturally. For a nice annotated 
collection of basic papers, see [4] . 

For applications we then need solutions of YBE. Many solutions have been found before, 
either by Lie-algebraic methods [5], or by using a specific ansatz [6]. In general YBE has 
iV 6 cubic equations for iV 4 unknowns, so even in the simplest, i.e. two-dimensional case, 
one has 64 equations for 16 unknowns. Such a set of equations is certainly too complicated 
for a brute force approach. In order to find the complete solution it is therefore necessary 
to simplify the problem by all means possible (without sacrifying generality). This is done 
using the inherent symmetries of the system. 



In writing out the two indiced object R we use the usual matrix notation. In two 
dimensions we have four 2x2 blocks, the second index pair denotes the location of the 
block: 



R 



/ pll p21 


pl2 p22 \ 

R u R u \ 




/ a b 


c 


d \ 


pll p21 

n 21 n 21 


pl2 p22 

n 21 n 21 




f 9 


h 


J 


pll p21 

n 12 n 12 


pl2 p22 

n 12 n 12 




k I 


m 


n 


\ pll p21 

\ rt 22 rt 22 


pl2 p22 , 
- fx 22 rt 22 / 




\ p q 


u 


v J 



(3) 



2 Symmetries 

The quantum YBE is invariant under the following continuous group of transformations 

R — k(Q ® Q)R{Q ® QY\ (4) 

where Q is a nonsingular N x N matrix and k a nonzero number. One representative is 
sufficient of all solutions related by (4), and during the solution process we will choose that 
representative which produces the simplest equations. We will use up this rotational freedom 
in a specific order, using the following parametrization: 
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In particular note the scalings with k and a diagonal Q. Let us define the scaling weight of 
a matrix element by w(R^) = k + / — i — j, then two nonzero elements of R with different 
weights can be scaled to one using A and k. 

The YBE has also discrete symmetries related to index changes: 

R ij -»■ R w ( 6 ) 

R% - R^jtl (indices mod N) (7) 
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In two dimensions and using the matrix notation (3), (6) corresponds to a reflection across 
the diagonal (= the usual transposition, also called the P reflection), (7) with n = 1 and 
followed by (6) corresponds to a reflection across the secondary diagonal (C reflection), and 
finally (8) corresponds to a reflection among the two central rows and among the two central 
columns (T reflection). 

3 Breakup into smaller parts 

Since the system under stydy is so large we must first split it by hand into several smaller sub- 
problems. [When we refer to specific equations, we use the following numbering: Ei 3+ 2i 2+ 4,i 1+ 8j 3 +i6j 2 +32ji-62 
R hJ2 Rl k\%R l k% ~ R tn R fik\R l k% The equations are written out explicitly in [2,3].] 

i) First we analyzed the equations by counting how many times each variable appears in 
them and it turned out that the corner elements d and p appeared most frequently. Thus it 
seemed to be good idea to transform so that d = 0. This can always be done: If p = we 
use the P reflection, else we use the B part in (5) to get 

d new := B 4 p + B 3 (f + k-q-u) + B 2 (a- g-h-l-m + v) (] 

+B(-b-c + j + n) + d, {) 



and since p^O we can always find a B so that d new = 0. 

ii) At this point we had d = and to keep it that way we take B = in subsequent 
transformations. When looking at the equations we found that one of them had a nice form 

E 22 := bc{f -k)+ jn{q - u) = 0, (10) 

and thus the problem would factorize into three parts, if we could transform so that f — k — 
or q — u = 0. This was accomplished as follows: 

If / = k already there is nothing to do, if / ^ k but q — u use C reflection to put f = k, 
and in both case take C = in Q (5). If both / ^ k and g^«we have after transforming 
with the C part 
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(q — u) new : = C 2 (b — c) + C(g — h + l — m) + q — u = 0. 



(f-k) new : =C 2 (j-n) + C(-g-h + l + m) + f-k = 0, 



We can now solve for C in one of the these (and use reflection (7) if necessary), except if 
j — n, b — c, h — I, g — m, f ^ k, q ^ u, which will become case C. 

iii) When we use d = 0, / = k in E22 we find that the problem splits into two big cases A: 
q = u and B: n = 0, q 7^ u. (The case j = can be T reflected (9) into n = 0.) When these 
assignments are used we find that some other equations factorize and the problem splits into 
six simpler cases (Al, A2, Bl, B2, B3, C). 

v) The two remaining continuous freedoms are related to scalings by A and k (4,5). If 
there are two nonzero elements with different weights they can both be scaled to unity. If 
we have elements of equal weight only one of them can still be scaled to unity. With the 
scalings we can split the problem further into a total of 33 subcases. 

4 Computer solution 

The best way to analyze sets of polynomial equations is by using Grobner bases. This is 
a systematic approach (Buchberger algorithm) to sets of equations and defines a canonical 
form in a given ordering of variables. It has been implemented in computer algebra systems. 
Since we just need solutions to a set of equations it is a good strategy to factorize the 
polynomials when possible and split the problem into smaller ones. Thus for each of the 33 
subsubcases we computed the factorized Grobner basis using the 'groebner'-package written 
by Melenk, Moller and Neun [7] for the REDUCE 3.4 [8]. The raw output contained repeats 
and subcases which were eliminated by a separate program. (In the newest version they are 
eliminated automatically). In the end we had 96 solutions to analyze. 

5 Canonical form 

Since many of the solutions obtained above can be transformed into each other it is important 
to bring them into a canonical form (using the continuous and discrete symmetries discussed 
before) for a final comparison. 

In order to define a reasonable canonical form let us consider the trace matrices of R: 
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Under (4) all of these transform according to 

r -> kQtQ- 1 . (13) 

The basic definition of the proposed canonical form is that the above trace matrices are in 
the Jordan canonical form. Since the trace matrices do not necessarily commute it is possible 
that they cannot be brought to the canonical form simultaneously, we will therefore work in 
the above order. If this requirement is not enough to fix the rotational freedom completely, 
we must look at individual matrix elements of R. (For a detailed algorithm for constructing 
the canonical form, see [1,2].) Using these ideas we were able to combine the 96 solutions 
into 23 cases using homogeneous parametrization [1]. 

6 Conclusions 

With the work [1,2] we know all solutions of constant quantum YBE in 2 dimensions. Most 
of them fit into the 8-vertex ansatz, but not all. One interesting observation is that all 
nonsingular solutions are either upper triangular, or have the property i?f- = unless i+j = 
k + l (mod 2). 

There are now many things that we can do with the results. For example, what kind of 
algebras do we get with the new solutions? And what kind of applications are related to 
them? 

As for extensions, next one could search for solutions with a spectral parameter. Pre- 
liminary studies indicate that this can also be done systematically and completely, but the 
work involved seems to be rather extensive. Another interesting problem is to go to higher 
dimensions. 



References 

[1] J. Hietarinta, Phys. Lett. A 165, 245 (1992). 

[2] J. Hietarinta, submitted to J. Math. Phys. 

[3] J. Hietarinta, in Proceedings of ISSAC'92, PS. Wang (ed.), p. 350, ACM Press (1992). 

[4] M. Jimbo (ed.), Yang-Baxter Equation in Integrable Systems, World Scientific (1990). 

[5] M. Jimbo, Commun. Math. Phys. 102, 537 (1987); V.V. Bazhanov, Commun. Math. 
Phys. 113, 471 (1987). 

[6] K. Sogo, M. Uchinami, Y. Akutsu and M. Wadati, Prog. Theor. Phys. 68, 508 (1982); L. 
Hlavaty, J. Phys. A: Math. Gen 20, 1661 (1987); ibid 25, L63 (1992); D.E. Radford, 'So- 
lutions to the Quantum Yang-Baxter Equation: the Two-Dimensional Upper Triangular 
Case' (Math. Dept, Univ. of Illinois at Chicago preprint, Jan. 14, 1991); S.-M. Fei, H.-Y. 
Guo and H. Shi, J. Phys. A: Math. Gen. 25, 2711 (1992). 

[7] H. Melenk, H.M. Moller and W. Neun, GROEBNER, A Package for Calculating Groebner 
Bases. (Included in the REDUCE 3.4 distribution package.) 

[8] A.C. Hearn, REDUCE User's Manual Version 3.4, RAND Publication CP78 (rev. 7/91). 



